Europhysics Letters 



PREPRINT 



Derealization and Heisenberg's uncertainty relation 

Gert-Ludwig Ingold(*), Andre Wobst, Christian Aulbach and Peter Hanggi 

Institut fur Physik, Universitdt Augsburg, Universitdtsstrafte 1, D-86135 Augsburg, 
Germany 



PACS. 05.60.Gg - Quantum transport. 

PACS. 71. 23. An - Theories and models; localized states. 



Abstract. - In the one- dimensional Anderson model the eigenstates are localized for arbitrar- 
ily small amounts of disorder. In contrast, the Harper model with its quasiperiodic potential 
shows a transition from extended to localized states. The difference between the two models 
becomes particularly apparent in phase space where Heisenberg's uncertainty relation imposes 
a finite resolution. Our analysis points to the relevance of the coupling between momentum 
eigenstates at weak potential strength for the derealization of a quantum particle. 



Introduction. - The derealization of two interacting quantum particles in a disordered 
potential has been the subject of intensive research since it was first addressed [|]-^). More 
recent work has generalized the Harper or Aubry- Andre model to study the behaviour of 
two interacting particles in a quasiperiodic potential |] |^j. Here, the interaction rather leads 
to a tendency towards localization. For finite densities, no clear indication of an interaction 
dependence of the phase transition in the Harper model has been found . 

Already in the absence of interaction, a quantum particle in one dimension exhibits a very 
different behaviour depending on the potential in which it is moving. For a periodic potential, 
the eigenstates are extended Bloch waves characterized by a quasimomentum. On the other 



hand, already a small amount of disorder suffices to localize the particle |10|. A different 
situation can arise when motion on a lattice is considered. Then, for a periodic potential 
incommensurate with the underlying lattice, a transition from delocalized to localized states 
occurs as the potential strength reaches a critical value. 

In order to gain more insight into the localization properties of a quantum particle in 
one dimension, we compare two lattice models, the Anderson model |l(]] and the Harper or 
Aubry- Andre model An analysis of the phase space properties of the energy eigenstates 

will reveal the dependence on the nature of the coupling between momentum eigenstates due 
to the potential in which the particle is moving. 

Random and quasiperiodic potential. - The Hamiltonian of the Anderson model is given 

by® 

H = -£>><» + 1| + \n+ l)(n|) +Wj2vn\n)(n\ (1) 
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where the Wannier states |n) are localized at the sites n = 1, . . . , L of a ring with periodic 
boundary conditions. The first term on the right-hand side describes the kinetic energy which 
defines the energy scale. The random potential of strength W is expressed by the second term. 
The distribution of on-site energies is determined by the coefficients v n distributed uniformly 
on the interval [—1/2; 1/2]. In the limit L — > oo the eigenstates of the Anderson model are 
known to localize for any nonvanishing potential strength W |Tc| ] . 

The second model of interest, the Harper model, is defined by the Hamiltonian 

H = J2 d n >( n + 1 l + \ fl + 1 )( n \) + A^cos(27r/3n) |n)(n| (2) 

n n 

where the random potential in ([!]) has been replaced by a quasiperiodic potential if the param- 
eter assumes an irrational number in the limit L — > oo. Then, the nature of the eigenstates 
depends on the value of the parameter A. For A < 2, all states are extended while they are 
localized for A > 2 || . 

For finite size systems, it is convenient to choose j3 = Fj_i/Fj where Fj_i and Fi are two 
successive Fibonacci numbers [jll|. In the limit of large systems /3 approaches the inverse of 
the golden mean, (\/5 — l)/2. With this choice of (5 the system contains L = Fi lattice sites 
with Fi—i periods of the potential. 

The Harper model possesses an interesting duality property || which becomes evident by 
transforming Wannier states \n) into new states 

|*) = L- 1 ' 2 ^ exp(i2irkl3n)\n) = iT 1/2 ^ exp fi2n^j^n^ \n) . (3) 

These are eigenstates of the momentum operator to eigenvalues fci^-i modi^. Neighbouring 
values of k therefore do not imply neighbouring momentum eigenvalues. 
With the transformation (|j) one obtains the dual Hamiltonian 



H= X - 
2 



(\k)(k + 1| + \k + l)(k\) + i J2 cos ( 2 ^ fc ) l fc ) W 



LA A k 



(4) 



By this transformation real and momentum space arc interchanged and the original potential 
strength A becomes inverted into 4/A. Comparison of (j|) and (|j) yields the self-dual point 
A = 2 which separates the parameter regimes of extended and localized states. In contrast 
to the nearest neighbour coupling in real space in (j|), the new Hamiltonian does not couple 
nearest neighbour momenta. The physical reason is that scattering by the incommensurate 
potential may change the momentum by a large amount. This coupling to quite different 
momentum values will be of central importance for our reasoning below. 

Real and momentum space. - After this discussion, the question arises why these two 
archetypical models behave so differently. In particular, what is the physical reason which 
allows a localization transition at finite value A = 2 in the Harper model? In order to answer 
this question, we first take a look at the structure of the wave function — ^2 n c n \n) 
expressed in terms of the Wannier states \n), which will provide information about the spatial 
extension of the state. An often used quantity is the inverse participation ratio in real space 

EMI 

^ = Ei c «i 4 - ( 5 ) 
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Fig. 1 - The inverse participation ratio in real space (full line) and momentum space (dashed line) is 
shown (a) for the Anderson model with L = 2048 and (b) for the Harper model with L — 10946. The 
curves represent averages over an ensemble of eigenstates of the respective Hamiltonian as explained 
in the text. 



Provided that J2 n l c «| 2 = 1> the inverse of this quantity indicates the number of lattice sites 
over which the wave function is distributed. A corresponding quantity 

ft = Ei d »i 4 ( fi ) 

n 

can be defined in momentum space, where 

d «= L ~ 1/2 E ex p( i2 ^£K ( ? ) 

The two quantities are depicted as full lines (real space) and dashed lines (momentum 
space) in fig. [l]a for the Anderson model of length L = 2048 and fig. 0b for the Harper model 
of length L = 10946. For the Anderson model, the curves represent an average over 50 disorder 
realizations with L/2 states around the band center each. For the Harper model, an average 
over all symmetric eigenstates has been taken ]l5| . 

For both models one observes a monotonously increasing inverse participation ratio in real 
space which corresponds to an increasing localization of the eigenfunctions as the potential 
strength is increased. Correspondingly, the inverse participation ratio in momentum space 
decreases with increasing potential strength, indicating a delocalization in momentum. The 
different limiting values of P x for strong potential reflect the fact that in the Anderson model 
the eigenfunctions localize at one site while in the Harper model two sites are occupied because 
we consider here symmetric eigenstates. 

While the overall picture is qualitatively the same for both models, we note an important 
difference which becomes apparent when the system size is changed. In the Anderson model, 
the transition from extended to localized states is smooth and notably shifts to lower values 
of W as L is increased. As a consequence, in the limit of infinite system size all states are 
localized if finite disorder is present. In contrast, for the Harper model one observes a sharp 
transition at A — 2 for sufficiently large system sizes. 
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Fig. 2 - The inverse participation ratio in phase space is shown (a) for the Anderson model with 
L — 2048 and (b) for the Harper model with L = 10946. The averages have been taken with respect 
to the same states as in fig. fil 



Phase space approach. The differences between the two models become much more 
distinct in phase space. This approach has the advantage of providing a consistent description 
for arbitrary potential strength Jl6| . The inverse participation ratio in phase space |Tij-]l9[| 

p = /^ir [g{x ' k)]2 (8) 

is based on the positive definite phase space density provided by the Husimi function |20[ or 
Q function pi| 

g(x ,k ) = \(x a ,k \t[))\ 2 . (9) 

Here, the state \tp) is projected onto a minimal uncertainty state centered around position xq 
and momentum fco- I n order to ensure equal resolution in the two directions of phase space 
we choose the width of the Gaussian as a = Ax = tJT/A^= l/2Ak. 

In fig. ^| we present the inverse participation ratio in phase space scaled with L 1 / 2 which 
is appropriate in the absence of a potential as well as for very strong potentials jl9| . For the 
Anderson model (fig. ^a), one obtains an increased inverse participation ratio at intermediate 
potential strengths implying that the eigenstates contract in phase space. As we will demon- 
strate below, the behaviour to the left of the peak is dominated by a contraction in real space 
corresponding to the increase of P x (cf. fig. |l|a) while to the right of the peak the decrease is 
dominated by the decrease of Pk- 

The qualitative agreement between the inverse participation ratios in real and momentum 
space for the Anderson and Harper model suggests that the same should hold true for the 
inverse participation ratio in phase space. This is even more so since the scenario just described 
for the Anderson model is consistent with the duality property of the Harper model where 
an inversion of the potential strength is accompanied by a transformation between real and 
momentum space. However, the results depicted in fig. ||b tell a different story. 

In contrast to the Anderson model, the inverse participation ratio P shown in fig. |^b ini- 
tially decreases with increasing potential strength up to A = 2. In this regime, the eigenstates 
therefore become more and more delocalized in phase space. Then, at A = 2, the phase 
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space distribution contracts and starts to become delocalized in momentum as A is increased 
further. Therefore, for almost all values of A, the phase space behaviour is dominated by 
the momentum component except for the transition which is dominated by the real space 
behaviour. 

The question now arises, why for weak potential the Anderson and Harper models behave 
so differently while for strong potentials they behave in the same way. The mechanism at work 
for weak potential can be considered to be responsible for the localization transition in the 
Harper model because the jump found in fig. ||b can only occur if the phase space distribution 
broadens as A is increased from zero. 

Uncertainty in phase space. Before addressing this question in detail, we recall the 
Gaussian smearing arising from the projection onto minimal uncertainty states which is in- 
evitable if a positive definite phase space density is required. As a consequence, the Husimi 
function only provides limited resolution which, as discussed above, we have chosen to be 
equal for the spatial and momentum components. The resolution is of the order of \[L sites 
which in the thermodynamic limit becomes small compared to the number L of sites in real as 
well as momentum space. It should be kept in mind, however, that for large L a phenomenon 
occurring on a fixed and finite number of sites cannot be resolved. These considerations would 
still hold, if we chose Ace and 1/Afc to scale according to L a with < a < 1. In contrast, an 
entirely different situation arises, if we keep the resolution fixed in one of the two phase space 
directions as is the case for the inverse participation ratios in real and momentum space, 
and (||). Then, no effects occurring in the other direction can be resolved even in the 
thermodynamic limit. 

Limit of strong potential. We are now in a position to answer the question raised 
above. It is useful to start by considering the limit of strong potential where the phase space 
behaviour of the Anderson and the Harper model agree. For infinitely strong potential, the 
eigenstates of the Anderson model are localized at one site while for the Harper model it is 
sufficient to restrict the discussion to one of the two sites occupied by a symmetric eigenstate. 
The kinetic energy may now be considered as a perturbation coupling in lowest order to the 
two neighbouring sites. The dominant contribution will be due to the neighbouring site closer 
in energy so that the problem reduces to the solution of a two-level system [pll . 

The energy eigenstates at very large but finite potential strength will be delocalized over 
two sites. While this will reduce the inverse participation ratio P x in real space, it will not 
affect the inverse participation ratio P in phase space which can only resolve spatial structures 
of size \f~L and larger. We recall, that this argument is independent of our particular choice 
of er since the absolute width of the minimal uncertainty state has to increase with increasing 
system size even though its relative width decreases to ensure a proper classical limit. 

The two original eigenstates which are coupled by means of the nearest-neighbour hopping 
term were localized on one site each and therefore totally delocalized in momentum space. 
Introducing the coupling, one finds a reduced spreading in momentum in order to ensure the 
orthogonality of the two states. This can easily be verified by considering the superpositions 
(|n)±|n+l))/-\/2 which in momentum space exhibit large scale density oscillations of period L. 
By means of the Fourier transformation the effect of the coupling on short real space distances, 
too small to be resolved in phase space, is turned into a large scale phenomenon. Therefore, 
in the momentum component the coupling is easily detected even with the finite resolution 
of the Husimi function. As a consequence, the behaviour of the inverse participation ratio in 
phase space is dominated by momentum and one finds an increase of the inverse participation 
ratio with decreasing potential strength as depicted in fig. |[ 
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Limit of weak potential. We now apply similar considerations to the regime of weak 
potential by starting from momentum eigenstates, i.e. states well localized in momentum 
analogous to the localized states in real space considered before. In the Anderson model 
the random potential will lead to a coupling among all momentum eigenstates. However, as 
before, the coupling between states close in energy, and therefore close in momentum, will be 
most effective. As a consequence, the role of position and momentum are interchanged and 
with the same arguments as above, we find an increase of the inverse participation ratio in 
phase space with increasing potential strength, albeit now due to the behaviour in real space. 
Within this perturbative treatment, we can readily understand the behaviour of P depicted 
in fig. |a. 

The situation is quite different for the Harper model where we may consider the dual 
model (^) for small A. The perturbation is now represented by the first term on the right- 
hand side of (Q) which couples to well-defined momentum eigenstates. However, as remarked 
below (|J), due to the scattering by the incommensurate potential these eigenstates in general 
do not correspond to nearest neighbour momenta. For most of the energy eigenstates, the 
momentum eigenstates to which the coupling occurs are far away on the scale of the resolution 
of the Husimi function. The resulting broadening of the momentum distribution leads to a 
reduction of the inverse participation ratio. In real space, on the other hand, the coupling 
will lead to short scale oscillations which are not resolved because of the finite resolution a in 
phase space. Therefore, the influence of the momentum component dominates and the inverse 
participation ratio decreases with increasing potential strength. 

Even in the Harper model there exist few particular states which couple to states close in 
momentum. Then, the inverse participation ratio will initially rise. However, the next order 
coupling leads to a distant momentum value and the inverse participation ratio in phase space 
will eventually decrease before reaching A = 2. 

It follows from this discussion that, in contrast to the Anderson model, the Harper model 
for both weak and strong potential is dominated by the momentum properties. It is only 
around A = 2 that real space becomes important. From the comparison of the Anderson 
and the Harper model we conclude, that the form of the coupling between the momentum 
eigenstates due to a weak potential plays a decisive role for the structure of the eigenstates in 
phase space and for the appearance of a delocalization-localization transition. 

As a further example we briefly comment on the Anderson model in two and three dimen- 
sions where the inverse participation ratio in phase space behaves very much like in the case 
of the Harper model (cf. fig. ^b) p9[ . In the marginal case of two dimensions, the tendency 
towards a transition is therefore clearly visible, even though the critical disorder strength 
vanishes in the thermodynamic limit |22] | and no true phase transition occurs. In contrast, in 
three dimensions the Anderson transition is recovered p3|, p4[. The main difference between 
the Anderson model in one dimension on the one hand and in two and three dimensions on the 
other hand lies again in the coupling of momentum eigenstates by a weak random potential. 
In the one-dimensional case, eigenstates close in energy are necessarily close in momentum. In 
contrast, in higher dimensions there may exist even energetically degenerate states far away 
in momentum, so that they can be resolved in phase space. 

Conclusions. - A crucial aspect of our discussion was the finite resolution available in 
phase space. This is in strong contrast to the ideal resolution available with the inverse 
participation ratio in real or momentum space, albeit only in one direction of phase space. 
As a consequence, there is no possibility to resolve the other direction even in the limit of 
large system size. Accepting Heisenberg's uncertainty relation and thus the finite resolution in 
phase space allows one to analyze the structure of the eigenstates in real as well as momentum 
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space and, as was demonstrated above, to obtain valuable information about the model of 
interest. 

This is corroborated by a recent observation by Varga et al. ]2j|, that instead of a full- 
fledged phase space calculation, one can alternatively make use of marginal distributions in 
real and momentum space. Even though the inverse participation ratios deduced from them 
resemble those defined in eqs. (|^) and (Q) a Gaussian smearing is again crucial. 

We therefore conclude that for the understanding of the localization properties of a quan- 
tum particle, where both position and momentum are relevant, the smearing in phase space 
called for by the uncertainty relation is not only necessary but also represents an essential 
ingredient of the physical argumentation. 

* * * 
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